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We present the calculation of K — > kk matrix elements of electroweak penguin operators i.e. Q-j 
and Qg. In the numerical simulation, we use the gauge configurations generated by the com- 
bination of Nf = 2 domain-wall fermions and DBW2 gauge action. From K — > n and K — > 
matrix elements on the lattice, we construct — ► 7T7T matrix elements at next-to-leading order in 
the chiral expansion by using recent analytic results. Renormalization factor of these matrix ele- 
ments are obtained by the non-perturbative renormalization technique in R1/M0M scheme. Brief 
discussion based on our results are made as well as a comparison with previous works. 
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1. Introduction 

In the efforts toward the theoretical treatment of the non-leptonic kaon decay, e'/e, the direct 
CP violation parameter of K — > nn, is one of the important focuses. This quantity is approximated 
to a linear combination of the K —> nn matrix elements of the local operators and Q% in the 
AS = 1 effective hamiltonian H^ S=1) = £?2i Wi(fl)Qi(fl) (where W t {n) is the Wilson coefficients 
obtained by perturbative calculations). However, calculations in lattice QCD (e.g. ^ Q]) resulted 
in unacceptable central values of e'/e due to a small magnitude of (nn ( /=0 ) |<26|^) compared to 
(nn ii-2) \Q%\K}. In these works, K — » nn matrix elements are obtained through the low energy 
constants (LECs) extracted from K — » n and K — > matrix elements on the lattice. Because this 
method is based on the lowest order (LO) chiral perturbation theory (ChPT), results could change 
significantly by taking the next-to-leading order (NLO) effects into account. 

For electroweak penguin operators Q ev/p = Qi$, there are limited number of NLO ChPT oper- 
ators due to their lower mass dimension than operators in other classes. By recent efforts ^ BJ BJ, 
there arose a prospection that the LO method can be extended to NLO for <2 eW p- In particular, the 
authors of Ref. [||] have shown that it is possible to extract all LECs needed to construct K — > nn 
at NLO from K — > n and K — > matrix elements in the framework of partially quenched ChPT 
(PqChPT). Following the line, we calculate K —>■ nn matrix elements of electroweak penguins on 
dynamical gauge configurations. 



2. Numerical simulation and analysis 



We use Nf = 2 gauge configurations generated by the RBC Collaboration ^ using domain- 
wall fermions with L s = 12 and M 5 = 1.8 and DBW2 gauge action with j8 = 0.80 in the 16 3 x 32 
box. There are ensembles with m sea = 0.02, 0.03 and 0.04, each of which 94 configurations are 
available, cr 1 = 1.69(5) GeV and am res = 0.00137(5) have been obtained in Ref. [@]. 

Electroweak penguin operators are defined as 



Q 7 = (s a Lfj,d a ) £ e g (qbR^qb), Q% = (s a L^d b ) £ e q (q b R^qa] 

q=u,d,s q=u,d,s 



(2.1) 



where L M = (1 -75)7^ and = (1 + 7>)7 M , and (e u ,e d ,e s ) = (2/3,-1/3,-1/3). Roman in- 

1/2 and 3/2 



dices represent color. Each operator is divided into two contributions with AI 

(1/2) (3/2) 

as <2ewp = <2ewp + <2ewp ■ We calculate correlation functions with the wall pseudo-scalar op- 
erator P(t) and the point operator <2ewp(x,f) for each m sea . With valence quark masses ra V ai = 
0.01,0.02,0.03,0.04,0.05, matrix elements are computed as a plateau in the ^-dependence of the 
ratios 
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(2.2) 



(2.3) 



(0|p(0Pt(o|0) 

where mp$ and mp$' are pseudo-scalar meson masses for degenerate and non-degenerate quark, 
respectively. Amp' PS ' is the amplitude of the PP-correlator. Locations of sources are (//,*/) = 
(5, 27) and we use t = 14-17 for the fit range. 
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One may write Q ewp as an expansion of ChPT operators with (8l, 8r) representation: 



Gewp = «88 (l^l ) + £c,-^ eWp) + 0(p% 



(2.4) 



r=l 



where ^ ewp ' are the NLO operators. Functions to which we fit our numerical results are ^ 
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(2.6) 



PS 



(2.7) 



where mj s 



m PS{m v 
„2 



=m sea )- For a set °^ non-degenerate valence quark masses such that m va n < 
m val2 , we define m z sS = m 2 s , (mvall=msca) and m\ s = '" P5 /( mval2 =m sca )- Pseudo-scalar decay constant f PS 
computed with degenerate quark masses is extrapolated to / in the chiral limit. After cancelling the 
1-loop divergences from the first term, there remain constant coefficients cj as LECs for NLO. Note 
that, in above functions, the number of LECs reduced from seven to five: {o^s, £"' 2 \ ^( 3 / 2 ) ) c r 4 , Cg} , 



where <g (1/2) = _ c r +£ .r +2 c r 3 + I0c r 4 + Sc r 5 and % ( 3 / 2 ) 



■ c r 2 + 4c r 4 + 4cj . One can easily check 



that these LECs successfully construct K ^ kk matrix elements at NLO 



3. Non-perturbative renormalization 

Before computing K —> Tin matrix elements, we renormalize lattice values non-perturbatively 
in RI/MOM scheme 0|. Numerical simulation for this step is carried out on 48,102 and 67 config- 
urations for m sea = 0.02,0.03 and 0.04, respectively. The main procedure is based on the method 
in Ref. [^]. We consider the renormalization factor of electroweak penguin Z as a 2 x 2 matrix to 
solve the mixing between Qj and <2s- The renormalization condition is 

z q 2 z ij r l *y uP2 ) = r%;, p.i) 

where T^ 11 is the amputated Green's function on the lattice in the momentum space. 

In the case of AI = 1/2, there is a contribution from the disconnected diagram in which an 
one loop contraction of <2; exchanges gluons with the spectator. For its reasonable implementation, 
we assign two momenta p\ and pj which satisfy \p\ \ = \p2\ = \p\ — Pi\ to each part of diagram 
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Figure 1: A i; - for AI = 1/2 (left) and 3/2 (right) in the chiral limit as a function of pjL. In each panel, 
diagonal (off-diagonal) elements are indicated by filled (open) symbols. 



so that all relevant momenta have same magnitude. In terms of lattice momentum pi att • 16/(27r) = 
(n x ,n y ,n z ,n t ), four sets (1,1,1, 1), (0,2,2,0), (1, 1,2,2), (2,2,2,2) and 4,8,8,4 equivalent partners 
to each are employed. We average results over the sets of momenta with same magnitude. 

Another issue for AI = 1/2 is the mixing with lower dimension operators. As discussed in 
Ref. this mixing with the two major operators can be subtracted as 

,(l/2)subt _ ^(1/2) , j , „ 2 



Q 



e! 1/2) + cj ■ (sd) + cf ■ s(D s + T>d)d 



(3.2) 



by determining the mixing coefficients c^ g and c 2 & . As a result of the procedure described in 
in Ref. we find the effect of subtraction is < 1% for the diagonal elements and < 6% for the 
off-diagonal ones. 



From (3A), the renormalization factor Z q 2 Z{j is the inverse of A ( y = Tr^g'T^ 66 ]. After the 
chiral extrapolation using the data with m sea = m va i, we show Ay in Figure [j] as a function of p^ att 
for AI = 1/2 (left panel) and 3/2 (right panel). Momentum dependences for pf M < 1.5 imply the 
contamination of low energy effects. To avoid this, we employ the data at the largest momentum 
(«Piatt) 2 = 2.467 for the rest of analysis. Given differences between two largest momentum are 
< 2a except for the 87-elements for both AI, we expect this contamination is under control. This 
small difference also moderates the worry that the 0(a 2 ) error is too large for the NPR technique 
to work correctly. 

After obtaining Z ; y by an alternative calculation of Z q , we move to MS with NDR through the 
perturbative procedure [§]. At H = a -1 , our preliminary results are 



z (M_=\/2) 
MS, NDR 



0.599(21) 
0.010(30) 



-0.059(38) 
0.557(93) 



7 (A/=3/2) 
"MS, NDR 



0.564(22) 
-0.121(21) 



-0.060(15) 
0.550(30) 



, (3-3) 



where only the statistical errors are considered. While the diagonal elements are w 40% smaller 
than the perturbative results similar results have been obtained by RI/MOM in Ref. [|l0|]. 

4. Results and discussion 



We fit the renormalized matrix elements to (2.5), (2.6) and u2/7J) simultaneously to ensure 
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Figure 2: Renormalized K — > tt matrix elements of 07 (left) and (right) as a function of m^ s [GeV 2 ]. In 
each panel, matrix elements with AI = 1/2 and 3/2 are presented for m sea = 0.02 (circles), 0.03 (squares) 
and 0.04 (diamonds). Fit curves are obtained from simultaneous fits to the PqChPT results. 



common values of Gi88 and c r 6 . This fit is repeated independently for each m sea . Using all data 
points (5 points for K — > % and 10 points for K — ► 0), we find reasonable quality of fit with 0.2 < 
# 2 /dof < 0.7. In Figure ^, we plot K — > % matrix elements of Qj (left panel) and Q% (right panel) 
and their fit curves for each m sea . 

By using LECs from each m sea and physical meson masses and decay constants, we compute 
K —> Tin matrix elements at NLO. For comparison, we also obtain the LO value by neglecting 
higher order terms with same a&g. Figure ^ summarises the results of K — > KK matrix elements in 
the physical unit with AI = 3/2 for Qj in left panel and Q% in right panel. The horizontal lines show 
the average over m sea 's with the weights of the errors. Results from previous quenched calculations 
Refs. [|l], ||, 1 1 , [12|] are plotted in filled diamonds for comparison. We observe large enhancements 
from LO to NLO sandwiching all of quoted results both for Q-] and Q$. 

Because NLO analysis of (7r7T( /=0 ) |<26|^) is not available currently, it is not realistic to esti- 
mate e'/e using our NLO result of (n7l^ =2 ) \Q%\K}. Instead of that, we could use the LO result 
from the NLO analysis in the estimation in the framework of LO. It should be noted that we obtain 
consistent values to our NLO results by the LO method employed in Refs. [[T], ||]. Therefore, as- 
suming similar gap we have seen between NLO and LO in the quenched calculations, there should 
have been a large enhancements from higher order effects in Refs. [|l], |2|. Then it is conceivable 
to obtain a fair value of e'/e by using (%7i/ I= 2) |<2s|^) a f ter removing large higher order effects 
by the NLO analysis, which is only possible for unquenched calculation. However, more careful 
study of the systematics must be done for actual conclusion. 
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Figure 3: Results of K — > nn matrix elements with AI = 3/2 for Qi (left) and Q% (right) in the physical 
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